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We review the magnetic form factor deduced by Delrieu from tiie Gorkov's equation for a Bardeen- 
Cooper-Schrieffer (BCS) type-II superconductor near its Bc2 phase boundary, i.e. when its magne- 
tization is small. A numerical study of the form factor, field map, and field distribution follows. The 
characteristics of the transition from the low-temperature BCS to the high-temperature Ginzburg- 
Landau vortex lattices is studied. The exotic shape of the component field distribution and the 
form factor at low temperature and as a function of the external field intensity are discussed. Our 
numerical work should be helpful for the analysing of small angle neutron scattering and muon spin 
rotation vortex-lattice data recorded for BCS superconductors and maybe other superconductors in 
the clean limit. 
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I. INTRODUCTION 

The bulk properties of a vortex lattice (VL) of a type II 
superconductor are studied experimentally, among other 
techniques, by magnetization, small angle neutron scat- 
tering (SANS), muon spin rotation (/zSR) and nuclear 
magnetic resonance (NMR) measurements. To extract 
physical information on the investigated compound mod- 
eling of the VL magnetic properties is required. This is 
usually done using either the London or the Ginzburg- 
Landau (GL) models.^- While the London model neglects 
the vortex cores altogether — which is acceptable for low 
fields only — the GL model accounts for them. Although 
the GL theory is strictly valid only near the supercon- 
ducting critical temperature at low field it turned out to 
be a good approximation for a number of classical BCS 
superconductors.^ The GL model is usually found to pro- 
vide a proper description of the VL thoughout the mixed 
phase for unconventional superconductors also, assuming 
the London penetration depth and the Ginzburg-Landau 
coherence length to be effective parameters.!^ This is fur- 
ther discussed by Landau and Keller .121 

One of the interesting phenomena predicted yet in 1972 
by Delrieu is the diffraction of the Gooper's pairs on the 
periodic potential induced by the VL."* He showed that 
for clean superconductors at low temperatures and fields 
close to the upper critical field i?c2 exotic behaviours of 
the VL may be observed due to the Gooper's pair diffrac- 
tion. The spatial field distribution around vortex cores 
obtains a conical shapes and the positions of the mini- 
mal and saddle point fields interexchange. As a conse- 
quence the probability field distribution shows a linear 
tail around the vortex core field. Nearly at the same 
time E. H. Brandt came to the same conclusion based on 
a nonlocal theory of superconductivity.'^ In the follow- 
ing publications he presented analytical and numerical 
results for the nonlocal VL behaviour at fields close to 
Bc2 with an arbitrary impurity scattering,!^ and later for 



a broad range of fields.'^ 

Writing a simplified Gorkov'P! integral equation for 
the Green's function in terms of a set of linear algebraic 
equations for the time and space Fourier components, 
U. Brandt et al. have been able to cornpute analyti- 
cally the density of states under high fields.*^ La ter o n the 
magnetizatiorfW] as well as the field distributiorP^ were 
also obtained by U. Brandt et al. and Delrieu, respec- 
tively, using the previously derived results for the Green's 
function. The clean and dirty limits were considered by 
Delrieu in his PhD thesis. 

It is worth to note that the Cooper's pair diffraction is 
a property of clean superconductors. Therefore, in most 
of high temperature superconductors these diffraction ef- 
fects may be observable. 

A VL field distribution which may exhibit a high- 
field linear tail was reported from NMR measurement 
on vanadium.!l2] Later on the existence of the tail was 
confirmed by Herlach et al. for niobium using fiSK 
measurements. They were performed at low tempera- 
ture and for an external field i?cxt relatively close to Bc2 ■ 
The linear tail is qualitatively different from GL model 
expectation for which a field cutoff should be present 

The NMR and the /iSR results seem to support Del- 
rieu's predictions, in particular the linear tail. SANS 
data may also be consistent with them.l^ However, after 
more than 40 years there is still no definitive experimen- 
tal observation of the predicted exotic VL at low tem- 
perature. Probably the difficulty of reading the Delrieu's 
works has prevented experimentalists to perform the re- 
quired combined SANS and /itSR measurements. Here we 
review this work and present a numerical analysis of the 
form factor, field map, and field distribution. 

The organization of this paper is as follows. Section 
recalls the geometry of a VL and introduces useful re- 
duced quantities. In Sec. Ill the physical principles be- 



hind the computations of the magnetization and form 
factor are given. The following sections, i.e. Sec. |IV|and 
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Sec. [Vj deal with the magnetization and form factor, re- 
spectively. In Sec. VI the component field distribution is 
discussed. The numerical analysis of the form factor, field 
map and field distribution is the subject of Sec. |VII We 



It is related to the lattice parameter of the VL: 



end up with a discussion and the conclusions in Sec. VIII 



Analytical and numerical details can be found in five ap- 
pendices. In particular, the last appendix shows that 
the form factor can be expressed in terms of a reduced 
number of parameters. This result allows in Sec 



VII 



to 

easily study the crossover from BCS to GL in the VL 
structures. 



II. GEOMETRY 

We assume Boxt to be applied along the Z axis of an 
orthogonal reference frame, with the vortex tubes of a 
type II superconductor in its mixed phase running along 
that axis. The VL is taken to be composed of equilateral 
triangles. Therefore in the direct space the unit cell is 
defined by the three vectors Vi = Xix, V2 = X2X -I- 
and V3 = z, where x and y are two mutually orthogonal 
unit vectors perpendicular to the Z axis, z is the unit 
vector of the Z axis, and {Xi,X2,l2} are coordinates 
with the relations X2 = Xi/2 and I2 = Since 
the vortex tubes are taken as straight, a VL point is 
labeled by the two dimensional vector 

R.p,9 = PVi + qw2 = {pXi + qX2)yi + qY2y. (1) 

Following Delrieu's convention a lattice point in the re- 
ciprocal lattice is specified by the reciprocal vector 

Km,?i = -h-v\ + = — [-/il2X + {mXi + hX2)Y] , 



with the VL unit cell area 

Sc = X1Y2 



(2) 
(3) 



The mean value of the induction is denoted B^{r) or 
in short, where r refers to a position in the direct 
space. We have introduced the magnetic fiux quantum 
$0 (4'o = 2.06783 X 10"^^ Tm^), and the two vectors 
and V2 which define the unit cell in the reciprocal lattice. 

As examples, with an obvious notation, we have 
Ri,o = Xi(l,0), Ro,i = Xi(l/2,^/2), and 
R_i,i = Xi(-l/2, V3/2). For the reciprocal lat- 
tice points, we compute Kio = [47r/(-\/3Ari)](0, 1), 
Ko,i - [4^/(y3Xi)](-y3/2,l/2), and K_i,i - 
[47r/(\/3Xi)](-vl5/2,-l/2). The complementary three 
points in the direct and reciprocal lattices respectively 
can be obtained by symmetry. 

We introduce for convenience two reduced quantities. 
We define the magnetic length 



Xi 




2.693 X A. 



(5) 



Unless i?cxt is close to the lower critical field Sci, the 
magnetization of a superconductor is negligible, and 
therefore ~ Bext- We shall express the form fac- 
tor with the unitless parameter nK„ ^ defined through 
the wave-vector scalar product 



TV- Ty'2 



From this definition, we derive the important formula: 



(6) 



'^K^.f. = i^^ + mh + h^) 



V3 



mh 



(7) 



III. PHYSICAL PRINCIPLES FOR THE 
DESCRIPTION OF THE MAGNETIZATION AND 
FORM FACTOR 

The computations of the magnetization and form fac- 
tor are based on an approximate form of the integral 
Gorkov's equation for the temperature Green's function 
Guji (r, r')P^ 

G„,(r,r') = GS,(r-r')- (8) 
I G^^(r-ri)l^(ri,r2)G!^„,(r2-ri)G„,(r2,r')d3rid3r2, 

where y(ri,r2) is the correlation function of the order 
parameter A(r): 



]/(ri,r2) = A(ri)A*(r2)exp A(l) • dl) 



(9) 



Note that F(ri,r2) describes the correlation function of 
the Cooper's pairs rather than the correlation of the elec- 
trons. A justification of the correlation function nature 
of V^(ri,r2) is given after Eq. 



13 We have introduced 



the Matsubara angular precession frequency ujf. 
uji = {2£+l)TTkBT/h. 



(10) 



The path of integration over the potential vector A is a 
straight line between ri and r2. Within the semiclassi- 
cal approximation the effect of the field on the material is 
entirely described by the phase integral in the correlation 
function and G^^{r—r') refers to the temperature Green's 
function of the normal metal and in the absence of a mag- 
netic field. The semiclassical approximation is valid if the 
spacing between the Landau levels is small compared to 
the sum of their thermal and collision broadenings, i.c^^^ 
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Here rnte is the level lifetime which accounts for the fi- 
nite electron mean-free path. It is futher discussed in 



Sec. VIII Neglecting the h/rnfc term in the ratio ex- 
pression, we compute T?.^ — 1 when T — 0.04 K for 
Scxt = 0.4 T. This Boxt value corresponds app roximately 
to Bc2 for a high quality niobium sample.^^ The theory 
for the form factor and resulting field map and field dis- 
tribution discussed here is therefore expected to be valid 
down to 0.04 K for simple superconductors such as nio- 
bium when the VL is composed of equilateral triangles. 

The correlation function has the periodicity of the 
VL with respect to the center of mass of a Cooper's 
pair, I.e. (ri + r2)/2.liSI Therefore it can be Fourier ex- 
panded. Neglecting the spatial variation of the induc- 
tion, a valid approximation in our case since we are in- 
terested in t he h igh-field VL for which the magnetization 
is negligible ,1^*^ 

^(ri,r2)= (12) 
X! ^K„,h(ri - ^2) exp [i'K„,,h ■ (ri + r2)/2] , 



m.h 



with 



^(r) = (-irVo(r-R 



(13) 



Since y(ri,r2) depends on the difference (ri — r2) and 
the center of mass of the Cooper's pairs, and not on 
ri and r2 individually, it is really a correlation func- 
tion. The Fourier components have a remarkable prop- 
erty: V^K„,JRm,/i) = (-l)'"^Vb(0). Therefore at the VL 
nodes, i.e. at position 'Rm,h, within a phase the cor- 
relation function of the order parameter has a common 
single value. In fact the phase changes from one posi- 
tion to the next nearest neighbors in a coherent fashion. 
Because of this long-range coherence, we expect Cooper's 
pair diffraction on the vortex cores. However, the diffrac- 
tion can only be partial since the coherent nature of the 
correlation is only active for a finite set of Cooper's pairs, 
i.e. the balistic ones with trajectory through the vortex 
cores. It is the extended correlation which is at the ori- 
gin of the slow power decay of the form factor at low 
temperature discussed latter on. This property does not 
hold near Bd since the spatial variation of the induction 
has been neglected in deriving Eq. [13) 

As the correlation function, the Green's function is pe- 
riodic with respect to (ri -I- r2)/2. In addition to the 
half sum of the coordinates, the two functions depend 
on the difference of the coordinates. To the difference 
corresponds the continuous conjugate Fourier vector p. 
Performing the 'K.rn,h and p Fourier transforms on the 
approximate Gorkov's equation given at Eq. [8) an alge- 
braic set of equations is derived. An approximate simple 
solution has been proposed in Ref. [S] based on the fact 
that the terms of the set with K,„ji ^ are negligible 
relative to the ones with J^m,h — . Latter on, it has 
been shown to be satisfactory only for the computation of 
thermodynamic quantities. It cannot be used to describe 
the dynamics. It is valid even at low temperature. ^'J^ 



The free energy of a superconductor can be written 
in terms of G^^^ (r, r').!^ Focusing again on the case for 
which -Bext is near 8^2, and therefore the spatial variation 
of the order parameter can be negl ected, and assuming 
the Abrikosov's vortex solution,SIini 



OO 



T 

Tr/2 

sin 6* ( — — hue ) d6 
V2a 



Ah\oJe\ 



We have defined the function 



a = a{9) 



A 



(14) 



(15) 



(16) 



hvp sin 6 

and the variable ue which is the root of the equation 

U( = 2hijJta + AQa'^iv{iue). (17) 

The Fermi surface has been assumed to be spherical. 
As we are using the BCS theory, we are working in the 
weak-coupling limit. The function v{z) is defined in Ap- 
pendix [X) 

In addition to the temperature, the free energy de- 
pends on four parameters: the critical temperature at 
low field Tco, the mean order parameter squared Ag = 
I A(r)p where the bar is for the spatial averaging, the den- 
sity of states at the Fermi energy per spin, volume and 
energy in the normal state Nq, and the Fermi velocity vp. 
The order parameter can be related to basic parameters 
of the superconductor. Assuming that the average value 
Ao is equal to the value of the order parameter in zero 
field, according to BCS 



Ao(0) ^ TT 

fceTcO exp (7) 



1.7639, 



(18) 



where Ao(0) is the value of Aq at T = and 7 the Euler- 
Mascheroni constant, i.e. 7 = 0.57722. So at low tem- 
perature only three material parameters are left if the 
spherical Fermi surface and weak-coupling aproximations 
are valid. 

In addition to Ao(0), other parameters characterize 
a superconductor. The Pippard-BCS coherence length 
^o(O) is related to An(0), 



eo(0) 



hvp 



(19) 



and in the clean limit - see for example Ref. [T] at page 
120, 



ecL = eo/0.96. 



(20) 



Here ^gl is the GL coherence length. We stress that 
Eq. [20] is derived in the T = limit. We note the GL 
relation 



Bc2{T) 



$0 



(21) 



It is also possible to derive information on the mean 
order parameter from the minimisation of the free energy 
written above. This leads to a formula of the Helfand- 
Werthamer typ^^^ which can be used to model Bc2{T) for 
an isotropic Fermi surface in the weak coupling approxi- 
mation. However, t hat ty pe of formula does not describe 
Bc2(T) for niobiumEsEl 

one of the metal for which the 
theory discussed here may apply. This is attributed to 
the strong anisotropy of the Fermi surface .I^^HU Hence, 
we shall not pursue any longer our discussion of Aq in 
terms of the free energy discussed in this work. 



IV. MAGNETIZATION 

While our main purpose in this report is the study of 
the form factor, here we discuss the magnetization M . 
This is justified since, thanks to Abrikosov, the relation 
between M and the induction is known for a temperature 
in the vicinity of Tco- Therefore the study of M gives us 
the possibility to check the validity of the formula for the 
induction, and therefore the form factor. 

We recall that M = -dF/dW. In terms of the 
free energy expression given at Eq. 15 dF/dB^ = 
— [a/{2B^)]{dF/da) and therefore we need to evaluate 
a{d/da){ui/a). This is done in Appendix [B| We finally 
obtain 



M = 



TTkBTNoAl °° 



2B^ 



E 



d{huJi) 



d9. 



(22) 

M can be evaluated numerically and compared to experi- 
mental data. Ourselves we shall used it in our discussion 
of the form factor. Interestingly, it can be drastically 
simplified in the T — limit. From Appendix [C] 



2BZ 



(23) 



This result will be used in Sec. |V] As expected, M is 
negative. 



V. FORM FACTOR 

Here we consider the field Fourier component B^ ^ 
which is usually called the form factor in the SANS lit- 
erature. It is related to the space-dependent induction 
component B^{r) by the Fourier relation: 



B^ir)= Sk„.. exp(iK„,,,. • r). 



(24) 



K„ 



with 



Bl^.. = - I S^(r)cxp(-zK„,,^ •r)d2r. (25) 



According to Refs. [4111] and correcting for misprints 
3^7VoA§ (-l)"'exp J 



MoT 



E 

1=0 



B^ 

tt/2 



sin(6i).g£(e')d6'. 



27rfcBrx 
(26) 



where we have introduced the auxiliary function 

iv {iui 4- mK„, h) + {iui — ifiK^ h) ^ 2i7j (iu^) dug 

d{hLUi 



1 



62 2° ■ 2 n V' (iUl) 



(27) 



It is convenient to rewrite the variable ug without the 



intermediate function a{9). Refering to Eqs. 16 and 17 
A 



A A^A^ 
ui = 2uje — + . „ ° . 2 J v{iue). 



(28) 



Therefore gi{0) is indeed a function of the angle 9 and 
the angular frequency uji. 

This formula for B^ ^ is extremely complicated. Re- 
markably, as M does, B^ ^ depends on three material 

^ 

parameters, i.e. Aq, Nq, and vp, and T and B . This 
suggests that simple relations may exist betweem them. 
They do exist in two limiting cases. Interestingly, Nq 
only appears as a proportionality parameter. 

In the proximity of T^q the asymptotic limit of the form 
factor is (see Sec. D 1 1: 



B 



\'mh 



exp 



(29) 



Recalling the definition of the Fourier transform of the 
induction, see Eq. [24] 

B^(r) = J2 ^i„..exp(iK,„,^-r).(30) 

This means that 
B^{r)=W~ (31) 
fio\M\ J2 (-l)™"exp(-n2,^^)exp(zK™,;..r), 

(m,/i)^(0,0) 

since M < 0, see Seclrvl This is the Abrikosov result.'^^l 



An exotic behaviour is found in the T — !■ limit. From 
Sec.lD2]for tt-k^.^ > 1, 



B 



K„ 



B^ 



In terms of the magnetization, see Eq. [23) it gives 



5k„. = -27rV^/io|M 



(32) 



(33) 
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Whereas ^ has a Gaussian wave-vector depen- 
dence near Tco , it decays much slowly with a wave- vector 
power law at low temperature. This is quantified in 
Sec. |VII| Hence, whereas only the low indices Bragg 
spots might be observed by small angle neutron scat- 
tering (SANS) near Tcq, more spots should be detected 
at low temperature. In fact, in addition to the usual six 
Ki Bragg reflections, Ki i and K2.0 spots have already 
been observed for niobium at 4.2 K for Bcxt parallel to 
the [111] crystal direction.'^ Only for that field direc- 
tion is the VL triangular, and hence is of interest here. 
Because of possible double Bragg scattering effects and 
large distortions in the VL, it is not easy to derive infor- 
mation from the published form-factor data on niobium 
to compare with the power law prediction. 



VI. FIELD DISTRIBUTION 

It is difficult to analyse analytically the spatial de- 
pendence of the induction. Here we shall focus on the 
field distribution as measured by /iSR and NMR tech- 
niques. We shall assume a disorder-free VL and the field 
width of the distribution to be small relative to B^ , i.e. 
only the distribution of the field component along Bcxt 
is measured.^ This latter condition can be checked to 
be fulfilled by looking at the computed distribution. For 
simplicity we shall write the distribution Dc{B^) without 
specifying that it does depend on i?oxt- 

Mathematically, the distribution can be expressed in 
terms of a two-dimensional Dirac function: 



DciB^)^J^ 5[B^{v)-B^] 



d^r 

Sc. 



(34) 



where the integral only extends over the unit cell. 
An important characterisation of a distribution is its 



variance: 



E 



(35) 



For the simple asymptotic form factor given at Eq. |32[ 
we derive for the standard deviation 



33/VoiVoA2 



B^ 



E 



1 



(36) 



5.756 



VII. NUMERICAL ANALYSIS OF THE FORM 
FACTOR, FIELD MAP AND FIELD 
DISTRIBUTION 



As shown in Appendix 



El 



the form factor B^ 



and c. The parameter a = —^j,qttN[)/^^2c/ B^ stands 
as a proportionality coefficient while b and c are dimen- 
sionless [5 — (A/tt^^)^ and c = A/^^ with the tem- 
perature and field dependent length scale parameters 
^f_ = hv-p / {'i.iTk^T) and = /i-wf/i'Ao, see Appendix 
E 1 . Interestingly, for Boxt — Bc2 & ^ 1 and c changes 
notably as T varies. As shown below, this enables us to 
easily discuss the physics of the form factor, field map 
and field distribution. 

Let us consider a classical BCS superconductor with 
Tco = 10 K. From Eq.[l8]we compute Ao(0) ~ 1.52 meV. 
Assuming Bc2(0) = 0.4 T, from the GL relation (Eq.[2T|) 
we find ^gl(O) = 28.7 nm 
eo^Ao(0)/;i 



Using the BCS relation 
and the approxima- 



19 



^'F — 4o7rAo^Uj/ft, see Eq. 

tion of ^0 in terms of ^gl given at Eq. [20) we esti- 
mate v-p ~ 2.00 X 10^ m/s. Let us assume for the 
second critical field the simple temperature dependence 
Bc2{T) — -Bc2(0)(l — t^), with the reduced temperature 
t = T/TcQ. The field and temperature dependencies of 
the gap are taken to be described by the compact formula 
Ao = Ao{T)VT^ = Ao(0)yr^vT^, with the re- 
duced field b — Bey^t/Bc2{T). Here we note that ^0 7^ 
since ^0 only depends on Ao(0) while is a function of 
Aq. Under these assumptions: 



/AAo 



Al{T)x{l-b) $0 Se2(T) 



1 eGi.iT) 1 - b 



1 - b 



0.110- 



(37) 



Above we used Eqs. T9]|2T The parameter c is field and 
temperature dependent: 



c ^ 1.1811 



(38) 



depends in fact only on three parameters noted a, b 



Interestingly, c is only written in terms of two reduced 
variables. Numerically, b = 0.110, 0.073, 0.047, 0.027, 
0.012, and 0.0011 for b = 0.5, 0.6, 0.7, 0.8, 0.9, and 0.99, 
respectively, c is linear in temperature in the low tem- 
perature limit and diverges as t — > 1. The temperature 
dependence of c for 6 = 0.5, 0.7, 0.9, and 0.99 is given in 
Fig.[T] In agreement with our discussion in Appendix |E 1[ 
while the c thermal dependence is quite strong and it 
gets large as the critical temperature is approached, b is 
weakly temperature dependent and has a negligibly small 
value in the Soxt Bc2 limit. 

The analysis of b and c suggests to start our discussion 
of the VL properties focusing on its c dependence near 
the & ^- limit. In Figs. [2] and [s] we show the contour 
plots of the spatial field distribution for b = 0.0011 {i.e. 
b = 0.99), and c = 0.01, 0.10, 0.30, 0.70, 0.90, 1.30, 3.00, 
and 20.00. Comparing the data at the top of Fig. [2] and 
bottom of Fig. [3) the field map, profile, field distribution, 
and form factor intensity are drastically different at low 
temperature and near Tcq. At c = Ccross — 0.9 a crossover 
occurs. This value of Ccross corresponds to T ~ 0.6Tco in 
agreement with a report of Brandt (this value depends 
on vf)^ 
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FIG. 1: Temperature dependence of the parameter c in re- 
duced temperature scale for four different values of the re- 
duced field h = Bcxt/Bc2, i.e. h = 0.5, 0.7, 0.9, and 0.99. 
The insert shows the same but on semilogarithmic scale in 
the vicinity of Tco. 



For c > Ccro ss we do observe behaviours expected in 
the GL regime.'^^'^S) a. minimum for the spatial field dis- 
tribution is found in between three neighboring vortices, 
while in between two neighboring vortices, i.e. at mid- 
point on the line connecting two neighboring vortices, a 
field saddle point is located. The signature of the mini- 
mum of the field is obvious in Dc{B^). The saddle point 
corresponds to the maximum in Dc{B^). As shown in 
the inserts of the field distributions, the intensity of the 
form factor is expected to be reduced for large indices. 

For c < Ccross the positions of the minimum and saddle 
point are reversed. The shape of the profiles are sub- 
stantially different. In the limit c — >■ the conical shape 
of the vortex profiles near the positions of the minimum 
and saddle point are particularly pronounced. This prop- 
erty is even more clearly seen at the minimum-field point 
as shown by the dashed line in the field profile. This 
feature reflects directly the weak power-law decay of the 
form factor which is a consequence of the Cooper's pair 
diffraction by the vortex cores. This can be qualitatively 
understood if we recall that the weight of the k harmonic 
for the infinite Fourier series of a triangle wave - a wave 
with a profile similar to the dashed line - is proportional 
to /{2k + 1)^. Note the alternating sign as in the 

expression in the T = limit, as well as the power 
law decay. As reflected by the different exponents of the 
power-law decays, this comparison is only qualitative. 
This is not surprising given the two dimensional nature 
of the field map. The triangle-wave model has only one 
dimension. At high field the BCS limit, i.e. T — >■ 0, and 
GL limit, i.e. T Tco, predict similar high-field Dc{B^) 
tails already for c > 0.70. But the low field features are 
still different. As noticed from Fig. |2j it is not required 
to go to extremely low temperature to observe Bragg's 
spots of large indices when Bext — >■ Bc2 . 

A linear high-field tail is predicted in Fig. [2] for 



Dc{B^). However, it is only expected at really low tem- 
perature, i.e. for c < 0.10, in contrast to observation.!^ 
The Dc{B^) shape near the minimum field is really dif- 
ferent near T = and Tcq. This difference remains to 
be seen experimentally. A signature of the crossover for 
Dc{B^) is in its sharp rise at low field without any shoul- 
der. Although never reported, it could be observed ex- 
perimentally. 

Up to now we have focused our attention on the physics 
very near the Bc2 phase boundary when the parameter 
b is negligible. Now we study it as we go out of the 
boundary. 

In Fig. [4] we present the maps near T = {c — 0.01), 
for c = 0.90, i.e. at the crossover temperature, as well 
as for c = 3.00. They are computed for b = 0.11, 0.44, 
and 0.99, corresponding to the reduced fields 0.5, 0.2, 
and 0.1, respectively. Focusing first on the maps at the 
top of the figure, i.e. for c = 0.01, we note that the BCS 
regime, i.e. when the Cooper's pair diffraction matters, 
disappears when leaving the Bc2 phase boundary. This 
is clearly seen as the saddle point moves in between two 
vortex cores, as expected in the GL regime. Physically 
the distance between the cores becomes so large that the 
Cooper's pair diffraction is no more operative. The recov- 
ery of the GL features of the VL appears at a lower field 
if the temperature is increased, as seen from the maps 
at the crossover temperature. As noted from the maps 
at the bottom, i.e. at high temperature, their properties 
are field independent. 

In Fig. [5] we show Dc{B^), the field profiles along the 
dashed and solid lines in Fig. |4j and the form factors. 
While at c = Ccross they are still weakly b dependent, this 
is no more the case when c = 3.00. This is obviously 
consistent with the field behaviour of the maps, as seen 
in Fig. |4] Confirming the strong field dependence of the 
maps at low temperature, i.e. c — 0.01, the shape of 
Dc{B^) and of the field profiles, as well as the amplitude 
of the renormalized form factors, are strongly field de- 
pendent. Only at high field is the behaviour in the BCS 
regime observed. 



VIII. DISCUSSION AND CONCLUSIONS 

As pointed out by U. Brandt et al. the approximate 
Gorkov's equation used above is valid in the region of 
fields near Bc2{T) where the magnetization vs field curve 
does not deviate appreciably from the straight \meW For 
superconductors with a GL parameter k ~ 1 this corre- 
sponds to a rather limited field range, however for k S> 1 
it covers a substantial part of the VL phase. The valid- 
ity of the model is related to two types of approximations 
which have been done for the derivations of the magne- 
tization and form factor. We shall focus our discussion 
on the latter quantity. The first type is inherent to the 
method and the second can be taken out if necessary. 

We first recall the two approximations of the first 
type. The form factor is computed with an approximate 
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FIG. 2: (Color online) (left panel) Contour plot M{x,y) = [B^{x,y) - BZ]/[B^^ - B^] for b = 0.0011 and four different c 
values, i.e. c = 0.01, 0.10, 0.30, and 0.70, from top to bottom, (middle panel) Field profile along the solid and the dashed 
lines shown in the left panel, (right panel) The corresponding component field distribution Dc(B^) is shown as a solid black 
line. For comparison, we present in dashed blue and dotted red lines Dc{B^) for c — >■ and c — > oo, respectively. In order to 
match the vortex core field B^^, the horizontal axis for Dc{B^) shown with the dotted red line (c — > oo) has been scaled {B^ is 
identical for all the curves). The insert shows the values of TZm.h ~ j^K^ ^/-Bki o| ^^^^ form factors in logarithm 

scale. To simplify the drawings, and without any lost of information, the dashed blue line for the Dc{B^) plots (corresponding 
to c — >■ 0) is presented only when c < 0.30. 
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FIG. 3: (Color online) The same as in Fig. [2] but for c — 0.90, 1.30, 3.00, and 20.00, from top to bottom, respectively. 



Green's function. First, the effect of the field on the 
function is only described with a phase integral. This is 
the widely used semiclassical approximation. The ratio 
defined at Eq. [Tl] has to be smaller than 1 for this ap- 
proximation to be valid. Second, the spatial variation of 
the induction is neglected. Therefore this cannot be valid 
if Bext is too close to B^i- 

We now discuss the second type of approximation. 



First, a spherical Fermi surface has been chosen. It 
should not be a problem to describe a superconductor 
with an anisotropic Fermi surface. However, it is proba- 
bly possible only numerically. Second, up to now the con- 
duction electron mean-free-path iaifp has been assumed 
to be infinite. Here we describe a method to account for 
the finite iy^ip value. 

In the case of an isotropic impurity diffusion, the ef- 
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FIG. 4: (Color online) Contour plots of the field distribution for three values of c (c = 0.01, 0.90, and 3.00), and for 6 
significantly larger than zero: h = 0.11, 0.44, and 0.99 (corresponding to reduced fields h = 0.5, 0.2, and 0.1, respectively). 



feet of these impurities can approximately be taken into 
account in the following way. First the Matsubara an- 
gular precession frequency tui in Eq. [TO] is substituted by 



^mfp/i'F- We identify 



Loi + l/(2Tinip) where 

introduced here with riifc used in Sec. |III[ Hence in the 
formula for the form factor lj^ has the meaning^'' 



uji = {21 + \)TikBT/h + l/(2ri„p). 



Secondly Aq has to be renormalised.'^ It is substituted 
by 



Aq 

1 - e{LOi) 



(40) 



with 



<^i)=n / (41) 



2 TinipWF Jo V ■'^F sin 9 



Note, lim^j^p^o £(w«) = 1- 



the scattering is not too strong, i.e. when Ti^p 



Since e{uje) > 0, when 
is suf- 
ficiently long, the renomalisation increases Ag. This 
means that the Pippard-BCS coherence length decreases, 
as expected. Consistent with the region of validity of the 
form factor expression given at Eqs. [26]and[27l the pro- 
posed renormalisation is only valid if i?oxt is not too far 
(39) from 3^2- In the expression the renormalisation 



occurs three times: twice explicitly and once through the 
variable Ui. As it has been done in the clean limit, it is 
possible to write the expression in terms of a re- 

duced number of parameters; see Appendix. |E 2| 



Eilenberger has derived approximate equations for the 
Green's functions for which a numerical method has been 
developed to solve them; see Ref. [28] and references 
therein. With this formalism the form factor can be com- 
puted with the effect of a finite ij^fp accounted for.^^ It 
would be worthwhile to reproduce the analytical results 
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FIG. 5: (Color online) Field dependence of the component field distributions, field profiles, and form factors at low temperature 
- c = 0.01 and b = 0.0011, 0.0275, 0.0733, 0.165, 0.44, and 0.99 (corresponding to 6 = 0.99, 0.8, 0.6, 0.4, 0.2, and 0.1, 
respectively) - for the left column, at the crossover temperature - c = 0.90 and b — 0.0011, 0.44, and 0.99 {b — 0.99, 0.2, and 
0.1, respectively) - in the middle column, and at high temperature - c = 3.00 and b — 0.0011, 0.11, and 0.99 (6 — 0.99, 0.5, 
and 0.1, respectively) - in the right column. 



described here with the numerical method and extend 
this study outside the region of validity of the analytical 
solution, i.e. far below Bc2{T). 

The linear Dc{B^) tail at low and high field and the 



for 



relative large amplitude ratio Tlm,h — 
Bragg's spots far from the center of the reciprocal space 
are consequences of the diffraction of the Cooper's pairs 
on the vortex cores. Hence, a priori they should also be 
observed even for non BCS s-wave superconductors. For 
these exotic properties to be found, measurements have 
to be performed on very clean single crystal supercon- 
ductors (which is usually the case for high temperature 
superconductors) at low temperature and for i?oxt suffi- 
ciently close to Bc2- 

For the diffraction of carriers on a periodic VL to mat- 
ter at low temperature £mfp should be substantially larger 



than the intervortex distance, i.e. d = A/^i„fp <Si 1. For 
d ^ 1 no diffraction takes place and the VL has the GL 
profile similar to that shown in Fig. [s] (see Appendix E 2 ) 
while for d — the superconductor is in clean limit which 
was discussed above. In Fig. [6] we show contour plots of 
the spatial field distributions, field profiles, component 
field distributions, and form factors for intermediate val- 
ues of J = 0.1 and 0.7 with fixed c = 0.001 and b = 0.0011 
(since 6 ~ and c ~ are optimal for the observation of 
the diffraction effects). For d < 0.1 the characteristics of 
the VL are similar to that of clean VL. With increasing 
d the exotic behaviour of VL gradually vanishes and at 
d ~ 0.7 the crossover takes place. This result differs only 
slightly from that obtained by E. H. Brandt using a non- 
local theory.l^ This crossover depends on the combination 
of &, c, and d. In Fig. [t] plots for the crossover condition 
are given as a function of h, c and d {b vs c for different 
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FIG. 6: (Color online) The same as in Fig.ji] but for d = 0.1 and 0.7 and fixed b = 0.0011 and c = 0.001. The different plots 
show clearly when the conditions for observation of a BCS-type of VL are achieved (i.e. for 6 ~ and c ~ 0) the transition 
from the BCS to the GL vortex lattice occurs at d ~ 0.7. For d — >■ the VL is in clean limit (see also Fig. |2| while for d — > oo 
it is of the GL type as in the high temperature range, i.e. c — )■ oo (see Fig. jsland Appendix E2l. 




FIG. 7: (Color online) Plots of fo vs c for d = 0.01, 0.1, 0.2, 
0.3, 0.4, 0.5, and 0.6 where the condition of equal fields at 
saddle and minimal points is satisfied (the crossover condition; 
see e.g. top panel of Fig. [3|. 



d values). This conclusive figure illustrates the following 
natural condition for the observation of the exotic VL be- 
haviour due to diffraction: the three length scales , f^, 
and £mfp should be significantly larger than the intervor- 
tex distance (here, ^'^ = hvp/{2TTkBT), — hvp/nAo, 
see Appendix IE 1|). 



In conclusion, we have reviewed the previous works 
of Delrieu on the exotic behaviour of the vortex lattice 
(VL) at high field and low temperature. It is the con- 
sequence of the Cooper's pair diffraction on the periodic 
VL potential. Analytical and numerical results for the 
magnetization and form factors are derived using the 
Green's function formalism. In agreement with previ- 
ous works of E. H. Brandt where a nonlocal theory of 
superconductivity was utilized (see e.g. Refs. |6I25| and 
references therein), we find a set of conditions for the 
observation of this exotic behaviour of the VL. Namely, 
the intervortex distance should be significantly smaller 
than each three length parameters: and fmfp (see 

the text). An expression for the standard deviation of 
the component field distribution has been derived. The 
results of a numerical study of the form factors {B^ ^ ) , 
field map and field distribution [Dc{B^)] have been pre- 
sented for a broad range of applied field i?oxt and cover- 
ing the whole range of temperatures from T = up to 
TcQ. In addition, effect of impurities was studied. This 
has enabled us to determine features which distinguish 
GL from low-temperature BCS vortex lattices. The be- 
haviours of the experimentally accessible Dc{B^) and 



quantities versus the normalized temper- 
ature and external field have been exposed. These results 
should at least apply to niobium for Bcxt || and 
maybe other classical BCS superconductors such as vana- 
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dium. This analysis will help in searching the exotic VL 
behaviour using the SANS, ^SR, and NMR techniques, 
since the Cooper's pair diffraction is not restricted to 
the BCS theory and the conditions of the diffraction pre- 
sented above can well be satisfied by most of the high 
temperature superconductors. 
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Appendix A: The v{z) function and the related 
Dawson's integral 

The magnetization and the form factor for the field 
near Bc2 are expressed in terms of the function v{z) which 
is related to the so-called complementary error function 
erfc(z)P 



v{z) 



exp(— r 

z - t 



-dt 



■ exp(— z^)erfc(— i^). 



We have the relation 

v{—z) — — 2V7i'j exp(— z^) — v{z). 



(Al) 



(A2) 



Here z and t are complex and real variables, respectively. 
In our case z = ix, where x is real, and in the asymptotic 
large x limit 



ivlix) 



1 



1 - 



1 



O 



1 



In general, we have the relation 

v'{z) ^ -2zv{z) +2, 



(A3) 



(A4) 



where v'{z) — Av{z)/Az. Combining the two previous 
equations, we derive 



and 



v' {ix) — l/x^, X 



iv"{ix) = -2/x^, X 



oo. 



(A5) 
(A6) 



The iv{ix) function is bounded as follows: 



^ (A7) 
We note the asymptotic limit of the Dawson's integral: 



r 1 

exp(-a;2) exp(t2)di = — , 



(A8) 



Appendix B: Evaluation of a{d / da){ui / a) 
Here we evaluate 



d fUi 



da \ a 



(Bl) 



which is required to derive the magnetization from the 
free energy. It is easily found that 



A 



dui ui 
da a 



(B2) 



To compute the dug/ da, we first note that according to 
Eq. [17] we can write 

— M£ + 2hu!ia + AQa^iu(w^) = /(it^ , a) — 0. (B3) 

This implies that 



da 



df/da ^huji + AQaiv{iue) 



df/due " 1 + Ala^v'{iui) 
where v'{z) — dv{z)/dz. This means that 
'U£v'{iui) + iv{iui) 



1 + Ala'^v'{iui) 
Now we note the relation 
d [uiiv{iui)] du£ 



(B4) 



(B5) 



d{hjji) d{hhJi) 



iv{iui) — U£v'{iui)] . (B6) 



Since 



due 



2a 



df/d{huj£) ^ 

dj/dui 1 + Ala^v'iiuiY 



(B7) 



d{hwi,) 
we derive 

d [uiiv{iui)\ 2a [iv{iui) — u^v' (iui)] 
dihtue) ^ 1 + Ala^v'iiue) ' 

Combining the previous equation with Eq. |B5[ we obtain 

Aq d [uiiv{iui)] 



(B8) 



A 



(B9) 



2 d{hLJi) ' 
Recalling the relation given at Eq. |A4[ we finally derive 



A 



^0 ■ ^ 9ui 



(BIO) 



Appendix C: Low temperature asymtotic limit of 
the magnetization 

We start from Eq. [22) When the temperature is very 
small, according to Eq. |10| it is justified to replace the sum 
2TrkBTJ2£ by the integral /^^ d{huj). Since, according to 
Eq. |A5[ v'{iue) vanishes when £ — > oo, 



M = 



7t/2 



sm{ey{iuo)de. 



(CI) 
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As a Matsubara frequency vanishes with the tempera- 
ture, 



Uq 



A2A2 



-iv{iuo). 



(C2) 



RecaUing Eq. 22 we find, as expected, that is pro- 



portional to the magnetization as written explictly at 
Eq.[29} 



Because we are focusing on the field region near Bc2, 
except for a small domain for which 9 can be small, we 
can take Uq = 0. Using Eq. A4 we then get v'{0) = 2, 



and finally derive the result written at Eq. [23| 

Appendix D: Asymptotic limits of the form factor 

Here we determine analytically two asymp toti c li mits 
of the form factor ^ starting from Eqs. 26 
We shall first study the high temperature limit. 



and 



27 



1. The behaviour near Tco 

When approaching Tco from below, B^ vanishes as 
does Bc2- Hence A is getting very large. Referring to 
Eq. [28j this means that Ui is large. This has two conse- 
quences. 

First, it is justified to consider the ug ^ ?^k„ ^ limit 
for the numerator of the fraction in ge{0). Recalling the 
Taylor expansion of a function. 



,2 



2iv (iui) 



(Dl) 



Secondly, let us now focus on the denominator, in partic- 
ular on the second term. Because ui is large, for iv{iui) 



we can use the first term of its expansion given by Eq. A3 
From Eq. [28] we then get 



U£ — 2uji 



A 



A2A2 



1 



vp sm ( 



Therefore, to a good approximation near Tco, 

A 



IbJe 



v-p sm ( 



(D2) 



(D3) 



Using Eq. 

AgA^ 
f?v\ sin^ 



5j we deduce 

v' (iue) = 



Aq 

2hujf 



{2£- 



^0 

l)27rfcBr 



(D4) 

Hence, since Aq vanishes on approaching Tco, the second 



term of the denominator in Eq. 27 becomes negligible 
relative to 1. 

Using the two previous results, we derive from Eqs. |26| 



and 



27 



the asymptotic behaviour of Sg- ^ near 



cO- 



B 



2TT^^ioNoAlkBT 



K„ 



E 



due 



BZ 

1x12 

sin{0)iv" (iui) 
d{hLo, 



(D5) 



2. Low temperature limit 

Since we are interested in this work by the limit for 
which the field is near Bc2, Aq is small. According to 
Appendix [A] v'{iue) is bounded. Therefore, except for 
small 9 values, we can neglect the second term in the 
denominator of the ge{9) function relative to one. As 
done for the study of the low temperature limit of the 
magnetization, we can substitute the sum 27rfcBT by 
the integral d{h(jj). This gives 



B 



^iVoMoAg (-1)™'' r/2 



B^ 



siii{9)he{0)d9, 



with 



(D6) 



hi{9) = exp 



(-"K.,„.) [/ 

J U 



iv{ix + in-K^ h)da^+ 



iv{ix — h)dx 



2 / iv{ix)dx 



(D7) 



In the first and second terms we use the new variables t = 
X + ?T.K,„ ^ and t = X — tik,,, ^ , respectively, and split the 
integration iv{it)dt into JJ*" iv{it)dt + iv{it)dt^ 
where Ua = uq + nj^^ ^ and Ua = uq — ^ ^'^^ the 
first and the second terms, respectively. The two inte- 
grals iv{it)dt cancel the third term. As a result we 
obtain the following relation: 



he{9) 



exp 

Uq 



iv{it)dt 



iv{it)dt 



(D8) 



Let us study the function he{9). Since the field is near 
Bc2 and we are at low temperature, uq ~ 0. Then setting 
Uq — and using Eq. |A2[ we get 



he{9) ^ exp (^-^n^^ ^'^ 



2V^ 



From Eq. \M\ 



exp(t2)dt - 2 



iv{it)dt 



(D9) 



iv{it)dt<-- ' [t- ^yt'^+4:/n)dt. 

(DIO) 

Hence the second term of hi{9), i.e. 
2 exp (^—nl^^ /q iv{it)dt, is negligible if ^ is 
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sufHciently large. Since the first term is proportional to 



the Dawson's integral, using Eq. A8 we finally derive the 
hg{9) asymptotic limit: 



hi{6) = %/7r/"K,„ h when nK„ 



(Dll) 



Combining this result with Eq. D6 we derive the 
asymptotic limit written at Eq. [32} 



Appendix E: The form factor in terms of a reduced 
number of parameters 



26fl28 



The original ^ expression given by Eqs 
depends on three material parameters Aq, Nq, vp and 
two experimental parameters T and i?oxt since ~ 
Bcxt- In the next subsection we show that in fact it is 
a function of only three independent parameters. Even 
more interesting, only one of these parameters has to be 
varied to study the region close to the Bc2{T) boundary. 
The second subsection gives a formula for B^ ^ when 
the electonic mean-free-path is a finite. The dirty limit 
is studied. 



1. The form factor in the clean limit 

It is easily shown that the formula for B^ ^ can be 
written as follows: 



B 



where 



h{0) 



(-l)"''exp(-n2, J - W2 



(El) 



l + h 



v' (iui ) 
sin^ e 



(E2) 

We have used the analytical expression of du(/d{huJi) 
written at Eq. |B7[ The m^'s are found to be the solution 
of the equation 



J2e+1) -ivliui) 



sm( 



sin^^ 



(E3) 



The proportionality coefficient a is in magnetic induction 
units. It is written as follows: 



a = — /xo7rA'oAQ27rfcBTz 



2A 



Mo^A^oAg||. (E4) 
B^ 



B^hvF 

We have also introduced the dimensionless parameters: 
6=(AM^f, (E5) 

and 

c = A/e, (E6) 



where we have defined the temperature and field depen- 
dent length scale parameters = Hvy / {"^Trk-QT) and 
^■^ = Hvp/ttAq, respectively. Therefore ^ has been 

written in terms of three parameters: a, 6, and c. Inter- 
estingly, b is vanishingly small when Aq — 0. Then the 



second term of the denominator of Eq. E2 is negligible 
and d becomes small. As a consequence, as expected, the 
form factor is also small. 

There is interest to study the temperature dependence 
of b and c near the phase boundary Bc2 (T) , in particular 
their asymptotes. We first note that Ag — ?► 0. In addi- 
tio n, by the definition of the second critical field given at 
21 A ~ ^GL since B^ ~ Bc2- Let us investigate the 
limit. According to Eq. 20 we derive A ~ ^o/0.96. 
Hence A is finite, and therefore 6^0. Since diverges, 
we also derive c — >■ 0. Concerning the T — >■ Tco limit, we 
note that A diverges as does ^gl- This first means that 
b is the ratio of two large numbers. _ Num erically we find 

Secondly, as 



Eq. 

T - 



VII 



6 <C 1. An example is given in Sec. 
is finite, c — > oo. 

According to this discussion, in the limit T ^ TcO and 
near Bc2, we can set c — c» and 5—^0. This means that 



U( = c : 



,(2^+1) 



sin 6 



(E7) 



and using Eqs. |Dl| and |A6| 

MO) - -n'^^^Jv"i^u,) ^ 2n^„,. g3(27+l)3 - 
Taking these results into account, we derive 



'"^"^ (E8) 



^K.,.=2|(-ir'^exp 



e=o 









1)-' 


X 


/ sin^ dd9 




Jo 



= 1.4024^(-l)™"exp(-n^^^). (E9) 

Here, we have used the results JJ^^^ sin'^ 0d6' — 2/3 and 
Y,^{2e+1)-^ = 7C(3)/8 = 1.0518, with ({s) denoting 
the Riemann Zeta function. As can be seen from the last 
equation, for T T^q the form factor converges to the 
GL solution. It is proportional to the factor d/c^. On 
the other hand, as it can be seen from Eq. |32j in the 
low temperature limit near i?ext — ^c2 c — >■ and 
& — >■ 0) B^ is proportionnal to d/c: 



B 



7nh 



(ElO) 



Thus, the two limits being proportionnal to a/c, it is 
convenient to use this field scale as units of field. 

In conclusion, considering the region very close to 
Bc2{T) and the limits near Tco and T = 0, we find that 
b is very small. On the other hand, c is large in the first 
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limit and negligible in the second. Consequently ^, 
and therefore the field map and distribution, is expected 
to strongly depend on c. This fact is used in Sec. |VII| 
for the study of the crossover from GL to BCS vortex 
structures. 



taken into account with only a single new parameter, 
i.e. d. 

In the dirty limit we have d 3> 1 . This implies 



c{2e+l)+d 



sm f 



(E16) 



2. The form factor for a finite mean-free-path Eg. |E14| converges to: 



In the case of a finite electronic mean-free-path, i.e. of 
a finite scattering rate l/rimp, the form factor depend on 
the four parameters a, b, c, and d as follows: 



where 



/((«) 



i=0 



(Ell) 



, 1 v'jiue) 
^ "(l-e«)2sin^( 



(E12) 



The Ui's are found to be the solution of the equation: 
c{2e + l) + d b 



Ul 



with 



(1 - £,)2 sin^ 



-iv{iui), (E13) 



'dr'\vU~^^^i±A±^]de. 



2 ./o \ 2sin6l 



(E14) 



Note, hm^_^^ ££ = 1. The parameters a, b and c are 
the same as in the clean limit case. We have defined the 
dimensionless scattering parameter. 



A 



1 €gl 



61/2 



(E15) 



which is approximately the ratio of intervortex distance 
to the electronic mean free path. Hence, relative to the 
clean limit case, the effect of impurities and defects is 



d 



2 7o 5(2^ + l) + d ci2l+l)+d 



Here we used Eq. 



A3 



and Jg^^^ sinxdx = 1. Therefore, 



21+1 + d/c 



\ 21 + 1 / 
The denominator of Eq. |E12| converges to: 



1 + 6- 



v' (iui) 



(1 - e,)2 sin2 . 



1 



and using Eqs. Dl and X6 



b 1 

52 (2^+1)2 



(E18) 



, (E19) 



[c{2i 



l)+d\ 1 



(2£+l) 



This means for the form factor: 
i?i„ ^^i(-ir''exp 

with 

4, 



(E20) 
(E21) 



3 



^ Y2l+\ + d/c)-\2l+l)- 
£=0 

(c2 + 6/(2£+l)2)-3 



(E22) 



Thus, as expected, no matter the temperature, if the 
scattering parameter d is large the VL properties are de- 
scribed by the GL model. 



^ M. Tinkham, Introduction to Superconductivity (McGraw- 
Hill, New York, 1996). 

^ I. Landau and H. Keller, Physica C: Superconduc- 
tivity 466, 131 (2007), ISSN 0921-4534, URL 
[htt p : / /www . scienced irect . com/science/article/ | 
I pii/S0921453407Q11100l 

J. E. Sonier, ReportsTrog. in Phys^70, 1717 (2007), URL 
|http : //iopscience . iop . org/6034-488'57fo7n/Rbl 

* J. M. Delrieu, J. Low Temp. Phys. 6, 197 (1972), 
URL http : / /www . sp ringerlink. com/ content/| 

P xvl2611085wj3vp6/?HUD=MP , 

° E. H. Brandt, Phys. Stat. Sol. B 64, 467 (1974), 



ISSN 1521-3951, URL |http : //dx . doi . org/ 10 . 1002/| 
pssb. 2220640207 
" E. H. Brandt, Phys. Stat. Sol. B 65, 469 (1974), URL 
http : //onlinelibrary . wiley . com/doi/ 10 . 1002/pssb . ~| 
[^2220650205/abstract 
^ E. H. Brandt, Phys. St at. Sol. B 77, 10 5 (1976), 
ISSN 1521-3951, URL [http : //dx . doi . or g/10 . 100271 
" pssb. 2220770109 ^ 
" L. P. Gorkov, Sov. Phys. JETP 36, 1364 (1959). 
^ U. Brandt, W. Pesch, and L. Tewordt, Z. Phys. 201, 209 
(1967), URL http://www.springerlink.coni/content/ 
I ?k=voiy.3ay.2820iy.29+p°j(3ay^209/i;^^ ' 



16 



^° U. Brandt, Phys. Lett A 29, 568 (1969), URL jhttp : //dx . 
r doi. org/10. 1016/0375-9601(69)90444-7 

^ J. M. Delrieu, Ph.D. thesis, University of Paris-Orsay 
(1974). 

A. R ung, Phys. Rev. Lett. 25, 1006 (1970), URL|http:| 
I // link.aps . org/doi/10 . 1103/PhysRevLett . 25 . 1006 
^■^ D. Herlach, G. Majer, J. Major, J. Rosenkranz, 
M. Schmolz, W. Schwarz, A. Seeger, W. Tempi, E. Brandt, 
U.Essmann, et al., Hyperfine Interactions 63, 41 (1990), 
URL http: //www. springerlink. com/content/?k=puby, 
r 3ay.28Hyperfiiie +interactioiis'/.29+vol°/.3a°/.2863°/.29+py. | 
r 3ay.284iy.29, 

A. Yaouanc, P. Dalmas de Reotier, and E. H. Brandt, 
Phys. Rev. B 55, 11107 (1997), URL http : //link. aps . , 
I org/doi/10 . 1103/PhysRevB ■ 55 ■ 11107| ' 
N. R. Werthamer, in Superconductivity, edited by R. D. 
Parks (Marcel Dekker, New York, 1969), vol. 1. 
G. Eilenberger, Z. Phys. 182, 427 (1965), ISSN 0939-7922, 
10.1007/BF01383120, URL http : //dx . doi . org/10 . 1007/ 
I BF 01383120. 

E. Helfand and N R. Werthamer, Phys. Rev. 147, 
288 (1966), URL |http : //link . aps. org/doi/10 . 1103/| 
PhysRev . 147 . 288 

S. J. Williamson and L. E. Valby, Phys. Rev. Lett. 24, 



1061 (1970), URL http : //link . aps . org/doi/10 . 1103/ 
PhysRevLett.24.1061, 



S. J. Williamson, Phys. Rev. B 2, 3545 (1970), URL http: 



I //link ■aps.org/doi/ 10.1103/ PhysRevB . 2 . 3545 

P. C. Hohenberg and N. R. Werthamer, Phys. Rev. 



153, 493 (1967), URL h ttp : //link. aps . org/doi/10 . 
P 110 3/PhysRev . 153.49 3 

K. Takanaka, Prog. Theor. Phys. 46, 357 (1971), URL 



|http : //ptp . ipap . jp/link?PTP/46/357 
A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 
J. Schelten, H. UUmaier, and W. Schmatz, Phys. Stat. Sol. 
B 48, 619 (1971), URL http: //onlinelibrary .wiley^ 
com/doi/10. 1002/pssb.2220480219/abstract, " ' 
"^^ A. Yaouanc and P. Dalmas de Reotier, Muon Spin Ro- 
tation, Relaxation, and Resonance: Applications to Con- 
densed Matter, International Series of Monographs on Ph- 
ysiscs 147 (Oxford University Press, Oxford, 2011). 
2^ E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995), URL 

http : //iopscience . iop . org/0034-4885/58/11/003 
^® J. F. Annett, Superconductivity, Superfluids and Conden- 
sates (Oxford University Press Inc., New York, Oxford, 
2004). 



^'^ G. Eilenberger, Phys. Rev. 153, 584 (1967), URL http: 
//link . aps . org/doi/10 . 1 103/PhysRev . 153 . 584 
P.Miranovic, M. Ichioka, and K. Machida, Phys. Rev. B 
70, 104510 (2004) , URL http://link.aps.org/doi/lo7l 
I 110 3/PhysRevB . 70 . 104510| 

P. Belova, M. Safonchik, K. B. Traito, and E. Lahderanta, 
Phys. Rev. B 83, 104518 (2011), URL http://link.aps7] 
r org/doi/10 . 1103/PhysRevB . 83 . 104518 

M. Abramowitz and I. A. Stegun, eds.. Handbook of Math- 
ematical Functions (Dover Publications, New York, 1970). 



